Degree of entanglement as a physically ill-posed problem: The case of entanglement 

with vacuum 
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We analyze an example of a photon in superposition of different modes, and ask what is the degree 
of their entanglement with vacuum. The problem turns out to be ill-posed since we do not know 
which representation of the àlgebra of canonical commutation relations (CCR) to choose for field 
quantization. Once we make a choice, we can solve the question of entanglement unambiguously. 
So the difnculty is not with mathematics, but with physics of the problem. In order to make the 
discussion explicit we analyze from this perspective a popular argument based on a photon leaving 
a beam splitter and interacting with two two-level àtoms. We first solve the problem algebraically 
in Heisenberg picture, without any assumption about the form of representation of CCR. Then we 
take the oo-representation and show in two ways that in two-mode states the modes are maximally 
entangled with vacuum, but single-mode states are not entangled. Next we repeat the analysis in 
terms of the representation of CCR taken from Berezin's book and show that two-mode states do not 
involve the mode-vacuum entanglement. Finally, we switch to a family of reducible representations 
of CCR recently investigated in the context of field quantization, and show that the entanglement 
with vacuum is present even for single-mode states. Still, the degree of entanglement is here difficult 
to estimate, mainly because there are N + 2 subsystems, with N unspecified and large. 

PACS numbers: 03.65.Ud, 03.65.Fd, 42.50.Dv, 42.50.Xa 



I. INTRODUCTION 

One of the fundamental postulates of quantum the- 
ory is that symmetry groups are represented unitarily. 
In principle, all unitary representations of groups are al- 
lowed. The choice of a concrete represenatation is dic- 
tated by the vàlues of 'quantum numbers'. Irreducibility 
of representations is not fundamental. Multiparticle Sys- 
tems are usually described by reducible representations 
since a tensor product of two irreducible representations 
is typically reducible. Restriction of physical analysis to 
an irreducible 'block' is not always justified, and super- 
positions of states belonging to different blocks may be 
acceptable. 

When it comes to quantum fields the main symme- 
try of interest is the one associated with the harmonic- 
oscillator Lie àlgebra of canonical commutation relations 
(CCR). The field is often assumed to consist of many, 
pcrhaps even infinitely many, oscillators. One often hears 
that the field is a system with infinitely many degrees of 
freedom. Another way of phrasing this is that the CCR 
corresponds to a system with infinitely many degrees of 
freedom, a fact implying by von Neumann's theorem that 
it has infinitely many inequivalent irreducible representa- 
tions. So even if we impose irreducibility as a constraint, 
which representation should we choose? Do all of them 
imply the same physics? Maybe different experimental 
configurations require different representations? More- 
over, if the field is indeed a multi-oscillator system, why 
should we impose the irreducibility constraint? - and so 
on. 

For the moment there is no clear answer to any of these 
qüestions, although an experiment that may distinguish 
at least between reducible and irreducible representations 



of CCR has been recently proposed [ij . The goal of the 
present paper is to show the freedom of choosing differ- 
ent representations of CCR is at the roots of yet another 
controversy. The problem of 'entanglement with vacuum' 
has led to the discussion whose never-ending character 
seems to be related to the fact the main question has 
never been precisely formulated. Its most recent mani- 
festation is the paper which will be the basis of our 
analysis. The physical configuration we discuss is the 
same as the one analyzed in 0]. 

Consider two two-level àtoms and a photon in super- 
position of two orthogonal modes. A Jaynes-Cummings 
type interaction couples each of these modes with a differ- 
ent atom. Initially the àtoms are in a product state, but 
finally they become entangled. This conclusion can be 
reached without any assumption about the form of rep- 
resentation of the CCR àlgebra that models the modes: 



(1) 



and I m = l\ n commute with everything (i.e. belong to 
the center of the àlgebra). Now, what can be said about 
the entanglement present in the photon state? Appar- 
ently, and this is the conclusion reached in a number 
of papers (cf. 0, 0, 01)» the state had to be entangled 
because the interaction is a product of two local opera- 
tions and thus cannot increase entanglement in the sys- 
tem. This is indeed true if one works with the irreducible 
representation constructed in terms of an infinite tensor 
product (we term it the 'oo-representation'). The oo- 
representation is very popular in quantum òptics commu- 
nity and is often treated as the representation of CCR. 
The problem is that this is only a representation and, 
according to the theorem of von Neumann 0, Ü| , there 
exists an infinite number of inequivalent irreducible rep- 
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resentations of this àlgebra. Another representation can 
be found in the Berezin textbook Q , so let us term it the 
l B- representation'. Now the conclusion from 0.0- EI is 
no longer true: The state of the photon is not entangled, 
but the interaction is no longer local, and this is why the 
final state is entangled. We first show how the two calcu- 
lations work in practice, and then perform a calculation 
in the reducible 'finite N representation' which, as shown 
for example in pj , reconstructs the Standard results with 
arbitrary precision for finite-time evolutions. Here the 
conclusion is again different: A single-mode state that 
is not entangled in the oo-representation becomes entan- 
gled with vacuum in the iV-representation if N > 1. 

So the conclusion will be that the very notion of en- 
tanglement is wcll defined only after having specified the 
mathematical representation of a state and, in particular, 
the number of subsystems it describes. But the charac- 
terization such as a 'two-mode state of a single photon' 



is not yet concrete enough to fix the mathematics of the 
problem. And this is the true source of the controversy. 



II. CALCULATION INDEPENDENT OF 
REPRESENTATION 

Consider the Hamiltonian 

H = iï f <g> A + R <g> = ^ f q \ . (2) 

A is an arbitrary operator, and R is the atòmic an- 
nihilation operator satisfying R 2 = 0, R\+) = \— ), 
R'\—) = |+), where |+) and |— ) are the atòmic excited 
and ground states, respectively. One can explicitly write 
the evolution operator without any need of specifying the 
form or algebraic properties of A: 



e -iHt = R ï R ^ cos ( t y/Ajtf} + RR t ^ cos (tVAiA) - itRÏ (g) sine (WAA^)A - itR® sine [WaTA)A^. (3) 

I 

Here sinca; = (sina;)/x and the square roots are well when applied to the Hamiltonian H = H\ + H2, 
defined since both AA^ and A^A are positive operators, 

no matter what A one seleets. Now take Ri = i?(g)l, R 2 = = RÏ ® i dk — Rk ® i «t , k — 1,2, (4) 

1 (8 R, and two annihilation operators a%, 02, satisfying 

the CCR àlgebra [a m ,al] = S mn I m . The formula Q implies e~ lHt = e -iHit e -iH 3 t^ 



e lHkt = R\Rk ® cos (ty a k al) + RkR\ ® cos (ty a\a k ) + tR[ (g) sine (ty a k al)a k - tR k <g> sine (ty ala k )a[. (5) 
Let |0) be a normalized vector satisfying a k \0) = 0, and let a = [a\ + a2)/v / 2· Then 

e -^*|-)|-) o t|0) = -^|-)|-)(cos(í^)aí|0)+cos(ív^)4|0)) 

(|+)|-) sine ( t yfh)h\Q) + \-)\+) sine {t^h)h\0)). (6) 
I 



This final state is in general entangled, although the de- 
gree of entanglement is difficult to estimate at such an 
abstract level. But for irreducible representations the 
central elements li, 1% are proportional to the identity /, 
and without loss of generality we can take I\ = 1% = I. 
Then 

e -ifft|_)|_ )a t| 0) = cosí|-)|-)at|0) (7) 

+^(|+)|-) + |-)|+))|0). 

The reduced density matrix describing the two àtoms will 
be important when we compare the predictions of irre- 



ducible representations with those of the reducible ones. 
It reads 

p = cos 2 í|-)|-)H(-| (8) 
+ Ísin 2 í(|+)|-) + |->|+>)((+|(-| + (-|(+|). 

These formulas are vàlid for all irreducible representa- 
tions of CCR. It is obvious that our state is maximally 
entangled at least at t — tt/2, but even now, as we shall 
see in the second of the following sections, we cannot 
conclude that the initial state a'|0) is entangled. 
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III. oo-REPRESENTATION CÀLCUL ATION 

The Hilbert space TLoo is an infinite tensor product 
of single-oscillator Hilbert spaces: Tíoo = The 
operators are 



(9) 



where / is the identity in Ti., and [a, w] = I- The vacuum 
is given by the infinite tensor product 



|0> = 10)®---® |0)® 



|0) c 



(10) 



Such a representation has all the mathematical patholo- 
gies typical of infinite tensor products. If the index 
n refers to a wave-vector in cavity, then a localized 
wavepacket must involve the entire infinite tensor prod- 
uct. However, it is typical to reduce discussions of two- 
mode problems involving a\ and 02 to 01 = a ® /, 
a-2 = I <S> a, and the vacuum to |0) = |0)|0). Then 
a\\0) = |1)|0), 4|0) = |0)|1), and 



ot|o) = i(|i)|o) + |o)|i; 



(11) 



This state is explicitly entangled so there is basically 
nothing lcft to prové. The statemcnt is made somewhat 
stronger by noting that the evolution operator in this 
representation is [! 



t(fi t ®o-_Rig)a t ) 



(12) 



We have ordered the Hilbert spaces in the tensor product 
so that the initial state can be writtcn as 



-L|-)|i)|-)|o) + -L|-)|o>| 



(13) 



-iHt 



is a product of two operators that act locally on the 
Hilbert spaces of 'the first atom plus the first oscillator', 



and 'the second atom plus the second oscillator'. In such 
a system the entanglement can be exchanged only locally, 
and since the final state at t = tt/2 involves nonlocal en- 
tanglement between the àtoms, the initial entanglement 
of the oscillators had to be the same. This is consistent 
with the obvious constatation that initially the oscilla- 
tors were entangled, and thus the state of light bchind 
a beam splitter involved maximal entanglement between 
modes and vacuum. 



IV. B-REPRESENTATION CÀLCUL ATION 



Consider the Hilbert space of column vectors of the 
form 



V 



(14) 



/ 



with the norm 



(*|*) 



E 

n=0 



dki...dk n \^ n (k 1 ,...,k n )\ 2 . (15) 



We apply the convention where all the degrees of free- 
dom of a photon are denoted by a single argument k, and 
dk is an appropriate measure. We assume the functions 
are symmetric in their arguments. The Hilbert space is 
therefore the direct sum of symmetrized tensor products, 
TÍb = X^o^® s ™- D en °te by P the projector on func- 
tions symmetric in the ^-variables, and let the functions 
/„ form an orthonormal basis in the Hilbert space Tí. 
The representation of CCR reads 



P 



/0 






fn(Pl) 





V2S( Pl - fci)/„(p 2 ) 



Pi 





ki)ó(p 2 



k 2 )fn(ps) 



(16) 



V 



a n and are related by Hermitian conjugation. It is 
understood that acting with a* n on |^) we perform inte- 
gration over the k variables. 



state is 



Let us now take two orthogonal functions, fi(p) and 
f 2 (p) , representing the wavepackets leaving a beam split- 
ter and localized at t = at its two sides. The vacuum 



|0) 



(17) 
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Acting on the vacuum with = (a\ + a\)/\/2 we find and 



° ,|0> " T2 





h(pi) + hipi) 




(18) 



-iHt\ 



->|-)«t|0) 



cosí I 



1 

Vi 



t \ 

+ h(pi) 




V 



siní 



V2 



/ 






w 



(19) 



/i(pi) + / 2 (pi) 




The initial state was 

v : y 

and the final one, at t = n/2, is 

í 1 \ 







eHA®H A ®H B (20) 



gHa^T-Ía^T-Íb, (21) 



where Tí^ is the Hilbert space of a two-level àtom. The 
dynamics is no longer a product one: The two atòmic 
Hilbert spaces get coupled through Hb, and there is no 
natural splitting of the single-photon subspace 7í oíTÍb 
into tensor products. The initial state is explicitly a prod- 
uct state of three subsystems. 

One can say that the same structure has the initial 
state in the oo-representation, but the third subsystem is 
in an entangled state. This is true, but the problem with 
the 73-representation is that the third subsystem does 
not have a natural decomposition into subsystems @- 
The distinction between the vacuum |0) and 'a vacuum 
|0) at one side of the beam-splitter' is undefmed in this 
representation since the vacuum is here unique, and 
corresponds to the vector (|17fl . An entanglement with a 
unique vacuum must be trivial. To the reader who is not 
yet convinced we propose an exercise: Try to prové that 
the initial state i|2U|) is entangled by cmploying any known 
entanglement measure. It turns out that the exercise is 
simply ill posed. 



V. TV- REPRESENTATION CÀLCUL ATION 

The idea of the representation is to make the num- 
ber of modes independent of the number of oscillators: 
Each oscillator is a wave-packet containing all the possi- 
ble modes, and the number TV of oscillators is a param- 
eter of the representation. In the weak limit TV — > oo 



the TV-representation reconstructs the results of irreduci- 
cle representations, but the formulas get automatically 
regularized in UV and IR regimes. Recently, in Q], it 
was shown that the representation can be, in principio, 
directly tested in cavity QED. The TV-representation of 
electromagnètic field operators was introduced in jlOl · 
and further analyzed and generalized in 0, 0, 0, Il4| . 

The representation is in the simplest case constructed 
as follows. Take an operator a satisfying [a, o'] = 1 and 
the kets | k) corresponding to standing waves in some cav- 
ity. We define 



a k = \k){k\ (g) a, 



(22) 



The operators satisfy [afc,a[,,] = Skk'Ik- Ik 1S m 
the center of the àlgebra. The fact that Ik is not pro- 
portional to the identity means that the representation 
is reducible. In our tcrminology this is the 'N = 1 
representation'. Its Hilbert space H is spanned by the 
kets \k,n) — \k)\n), where a^a\n) = n\n). A vacuum 
of this representation is given by any state annihilatcd 
by all etfc. The vacuum state is not unique and belongs 
to the subspace spanned by |fe,0). In our notation a 
N = 1 vacuum state reads |0) = J2k Ok\k, 0) and is 
normalized by J2k 1^1 2 = Sfc = 1. The probabil- 
ities Zk = \Ok\ 2 play an essential role in the formal- 
ism, and of particular importance is their maximal value 
Z = max/jjZ/j}. In relativistic theory Z is a Poincaré 
invariant that additionally parametrizes the representa- 
tion and plays a role of renormalization constant. For 
N > 1 the representation space is given by the tensor 
power 7£ = TL® N , i.e. we take the Hilbert space of N 
(bosonic) harmònic oscillators. For large N the effective 
paramcter that controls the representation is the product 
NZ. Lct A : H -> H be any operator for N = 1. We 
denote = I®^ 1 ) ® A ® 7®^-»), A^ : H -> H, 
for 1 < n < N. For arbitrary TV the representation is 
dcfincd by 



[Ofc>oI'] 



^E4"\ L^E^, (23) 

(24) 



n=l 
-k ■ 



n=l 
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The iV-oscillator vacuum is the TV-fold tensor power of The spectral projectors E n (s) are in the center of CCR. 

the N = 1 case, a kind of Bose-Einstein condensate con- Defining a n (s) = a n E n (s), 
sisting of N wavepackets: 

|0) = |0> ® • • • <g> |0) = |0)® w . 



(25) 



Let us denote çtj = a fe , a 2 = a^. . In order to apply the 
formula © we first spectrally decompose the operators 

L 1 =L kl ,L 2 =L· 2 ÏÏM, 



a n (s),a m (s ) T ] = (s/N)S mn S as 'E m (s), (27) 



N N 

^ = S l^ s )' ^ = l^ 5 )- ( 26 ) we rewrite ® as 



s=0 



J 




As before a — (ç^ + a 2 )/V^, and Aí is chosen to make 
Wa^O) normalized: 

|^| 2 (0|aat|0) = \M\ 2 (Z 1 + Z 2 )/2 = 1, (29) 
(OlOfcoilQ) = (0|I fe |0) = (0|/ fc |0) = Z k . (30) 

The formulas look much more complicated than those 



from the irreducible representations, but they turn out 
to imply the same physics if one appropriately defines the 
limit iV — * 00. The limit cannot be performed directly 
(i.e. in the strong sense) at the level of (|28|l . It has to 
be computed in a weak sense, for example at the level 
of averages or the atòmic reduced density matrix, which 
reads 



1 N I 

b;Ei->i-><-i<-i- 2 (V^) 



Z! + Z 2 



s=0 
N 



S 

N 



N 



N 



Z s 1 (l-Z 1 ) 



N-s 



s.s'—O 

Esin (t\ / ^ \ / — sin ( 
^ ,_ V V NZJ V N \ 



8,S'=0 

N 



1^X™\ 1 a!nz)n- 



N 



Z S 2 {1 - Z 2 ) 



N-s 



N 



Z s 2 {l-Z 2 ) 



N-s 



N 



, \Z{Z S 2 {1-Z 1 -Z 2 ) 



N-s-s' 



N 



, \Z{Z S % {1-Z 1 -Z 2 ) 



N-s-s' 



(31) 



6 



The law of large numbers for multinomial distributions implies 
Poc = lim p 

Z 2 o f . I Z<l 

■ cos 



^Tz~ 2 C0S \1Y 



Zi + Z 2 



Zi + z 2 
-)|-)(+K- 



Zi 



SUI t\ , , 

Zi + z 2 V V z ) V Zi + z 2 



__ sin (ty^)|- 



Z ' 2 M1, ('\ §)< H 



■M+X-K+l 
f |-)l+)(+K-|) 



(32) 



r 



A 2 

1 - Probability distribution Z k =10*1 



Z, Z 2 




Infrared #1 *2 Ultraviolet 



Quantum optical regime 

FIG. 1: Qualitative properties of the the vacuum probability 
distribution corresponding to the vacuum wave function in the 
iV-representation. For wave vectors in the 'quantum òptics 
regime' we expect Zk/Z ~ 1, where Z is the maximal value 
oïZ k . 



Various arguments based on comparison of fields quan- 
tized in iV- representations with those quantized by means 
of irreducible representations suggest that the vacuum 
probability distribution Zk is qualitatively of the form 
shown in Fig. 1. If the two wave vectors fci, k 2 cor- 
responding to a± and a 2 are assumed to belong to the 
plateau region (i.e. Zi = Z 2 = Z) then becomes 
identical to (JSJ). The dynamics looks as if it were reg- 
ularized by the cut-off function \k = \fZkj~Z if Z^ are 
in the IR or UV regimes. Let us also note that the dy- 
namics is not generated by H but by its renormalized 
version H / ^/~Z, a property equivalent in quantum òp- 
tics to bare charge renormalization. Internal consistency 
of calculations based on JV-representations requires that 
N < 00. Simultaneously, as discussed in detail in Q], for 
any finite-time evolution one can choose such an N that 
the evolution will be indistinguishable, within some given 
error bars, from the result based on irreducible represen- 
tations. Taking t — 7r/2 one can always choose a finite 
N guaranteeing that (|31f) is arbitrarily close to the maxi- 
mally entangled atòmic state © . Paradoxically, we have 
this agreement for a large N, and not for N — 2 which 
one might naively expect on the basis of the experiences 
with the oo-representation. 

Concerning the entanglememnt with vacuum of the ini- 



tial state, the iV-representation reads explicitly 
at|0) = -L(at|0)|0)...|0) + ··· + |0)...|0)at|0)). 

Here a = (ai +a 2 )/^/2 with eti and a 2 defined for N — 1. 
This state is highly entangled, but it would be entangled 
also for a single-mode problem, since 

fltlfl) = -L(at|0)|0)...|0) + ··· + |0)...|0)al|0)) 

even though its oo-representation is a product state. 
a{ 1 0) is similar to the oo-representation state 

* (|l)|0)...|0> + ··· + |0)...|0>|l)) 

representing a superposition of N different modes, but is 
'more entangled' since each N = 1 component a[\0) is a 
superposition of all the wave-vector degrees of freedom. 

VI. CONCLUSIONS 

It is not clear which representation of CCR should be 
used to model quantum electromagnètic fields. The von 
Neumann theorem guarantees that there exists an infi- 
nite number of different irreducible representations, and 
we have not found a single quantum òptics problem that 
required the use of the oo-representation and not, say, 
the Berezin one. Even the most spectacular problems 
involving entanglement, such as teleportation or viola- 
tion of Bell's inequalities, turned out to be solvable at 
a representation independent level. To the best of our 
knowledge all the formulas given in the recent reviews 
on linear òptics quantum computation [Ï5l IÏ6| can be 
translated into a representation independent formalism, 
and then into the i?-representation, say. Accordingly, it 
is absolutely irrelevant which irreducible representation 
we select for our explicit calculations. Thcsc remarks 
apply, in particular, to all the papers on entanglement 
with vacuum and nonlocality of a single photon where 
the calculations were performed in the oo-representation 

ElHHÜEïlillülHlil ( see however 0)- 



7 



No single experimental prediction would be affected by 
switching to the B-representation, but the discussion of 
entanglement would be changed, of course. 

The paper j2Ïj whose intention was to clarify the status 
of entanglement with vacuum has increased the confusion 
by making a distinction between 'entanglement in config- 
uration space' and 'entanglement in Fock space', but the 
Fock space was identified with the oo-representation. It is 
evident that any irreducible representation of CCR leads 
to the notion of a Fock space, i.e. a direct sum of Hilbcrt 
spaces constructed from a vacuum state by means of cre- 
ation operators. The conclusions of [27( would not be 
possible to formulate if the authors discussed the prob- 
lem in the B-representation of the Fock space. 

The situation changes if we consider the reducible rep- 
resentations. As opposed to irreducible representations 
here even a single-mode state of light is entangled with 



vacuum if N > 1. The A^-representation is in principle 
experimentally testable If we managed to determine 
the physical representation (representations?) of CCR, 
the discussion of entanglement in optical systems might 
finally become physically wcll posed. 
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